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I. INTRODUCTION 


The possibility of discovering the so called quark-gluon 
plasma at RHIC (Relativistic Heavy Ion Collider) or in 
the future LHC (Large Hadron Collider) has triggered 
a lot of effort in the theoretical description of ultra- 
relativistic heavy ion collisions. The modern view of 
these collisions was set mainly by Bjorken Q| and it is 
based in the hydrodynamic model already considered by 
Landau j|] . In the last stages of the evolution of the cen¬ 
tral rapidity region the hadronic fluid is made mainly of 
pions. It has been argued that this expanding pion gas 
can reach thermal equilibrium much more efficiently than 
chemical equilibrium since at low energies pion interac¬ 
tions are mostly elastic . In order to have an appropri¬ 
ate phenomenological description of the expanding pion 
gas it is needed to set the proper hydrodynamic equa¬ 
tions. In the pioneering work by Bjorken the hadronic 
gas (both in the confined and in the quark-gluon plasma 
phases) was assumed to be a perfect relativistic fluid 
and accordingly viscosity and heat conductivity were ne¬ 
glected. In order to check the validity of this assump¬ 
tion in the different regimes of the pion gas expansion 
it is needed to compute the corresponding transport co¬ 
efficients as a function of the temperature and density. 
This could make possible to know in what physical situ¬ 
ations they are relevant and the ideal fluid equations are 
not appropriate any more. In addition the expressions 
obtained for these coefficients can be used to set the cor¬ 
rect relativistic version of the Navier-Stokes equations 
which applies when dissipative effects must be taken into 
account in the pion fluid dynamics. Curiously enough, 
the computation of these coefficients is also interesting 
from the fundamental point of view, since it can be done 
completely from first principles. 
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In this work we will illustrate this by computing the 
viscosity of the pion gas at low energy and density. Our 
computation will rely just in standard kinetic theory and 
the chiral symmetry properties of the strong interaction 
and thus our results can be considered model indepen¬ 
dent. In fact our results will depend only on the pion 
mass Mtt and the pion decay constant A/. The main 
points on which our computation is based are the follow¬ 
ing. First of all our computation is going to apply only 
at low temperature and density. At low temperatures 
hadronic matter is confined into hadrons and at very low 
temperatures the only excited modes are pions for van¬ 
ishing barionic number ng density. Notice that ub ~ 0 is 
supposed to be the case in the central rapidity region of 
heavy ion collisions. Moreover the additional low-density 
condition allows us to consider the non barionic hadronic 
matter as a pion gas. 

At low temperatures most of the pion interactions oc¬ 
cur at low energies. The low energy pion interactions are 
completely determined by the chiral symmetry pattern 
of the strong interactions. In the chiral limit (vanishing 
u and d quark masses) the strong interaction dynam¬ 
ics is invariant under the SU{2)b x SU{2)b chiral trans¬ 
formations. However this symmetry is supposed to be 
spontaneously broken to the diagonal or isospin group 
SU{2)l+r, being the pions the corresponding Goldstone 
bosons. In the real world quarks masses are different 
from zero and then this symmetry breaking pattern is 
not exact resulting in a small value (when compared with 
other hadrons) for the pion mass M^. In any case chiral 
symmetry determines completely the low energy pion in¬ 
teractions in terms of and (Weinberg low energy 
theorems). The corrections appearing at higher energies 
can be computed in a wstematic way by means of Chiral 
Perturbation Theory M . In the chiral limit the pion scat¬ 
tering amplitude goes to zero in the low energy limit. Out 
of the chiral limit it is proportional to M^/F^. In any 
case the interaction is small in this regime. Therefore we 
arrive to the conclusion that non barionic hadronic mat¬ 
ter at low temperature and density can be described as 
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a weakly interacting pion gas. Another important obser¬ 
vation is that at low temperatures (energies) most of the 
pion collisions are elastic. This fact implies in particular 
that the number of pions is conserved. Thus it is possible 
to introduce the corresponding pion chemical potential /i 
(not to be confused with the more usual chemical poten¬ 
tial associated to the barionic number). This makes it 
possible to consider the pion gas at thermal equilibrium 
at temperature T for different values of the chemical po¬ 
tential or what it is the same, different pion densities, 
and not only the case of chemical equilibrium correspond¬ 
ing to /i = 0 . 

According to the previous discussion the simplest state 
for non barionic hadronic matter at low temperatures 
and densities is just a gas of weakly interacting pions in 
thermal equilibrium which could be described in a first 
approximation with a Bose-Einstein distribution. In ad¬ 
dition, at very low temperatures the average pion energy 
is much smaller than the pion mass allowing for to have a 
non relativistic description of the pion gas. As discussed 
later, the low density condition will prevent the forma¬ 
tion of the Bose-Einstein condensate out of the pions, 
thus making the description of the hadronic matter in 
the mentioned regime especially simple. 

More generally it is possible to consider the case when 
the equilibrium is only local. In this case the tempera¬ 
ture, density (or pion chemical potential) and the overall 
velocity are space-time functions. In this situation it is 
possible to have an hydrodynamic description of the sys¬ 
tem as that considered in the Bjorken picture of heavy 
ion collisions. A departure from local equilibrium will 
give rise to dissipative processes like viscosity, which is 
the main topic of this work. The evolution of the dis¬ 
tribution functions outside equilibrium can be studied 
by using the Uehling-Uhlenbeck equation, which is the 
quantum version of the Boltzmann transport equation. 
The computation of the transport coefficients such as 
viscosity or heat conductivity requires the solution of the 
transport equation for different specific perturbations of 
the local equilibrium distribution. 

As it was mentioned above, in this work we are inter¬ 
ested in the computation of the viscosity of the hadronic 
fluid with the approximations discussed in the previous 
paragraphs. In fact it is possible to have a good estima¬ 
tion of this magnitude by neglecting quantum effects. In 
this case the equilibrium distribution function is just the 
Boltzmann distribution instead the Bose-Einstein distri¬ 
bution. Moreover the elastic scattering pion cross section 
is constant at low energies and as we will see it is given 
by 


23 Ml 
384 ttFI 


( 1 ) 


In addition, at low temperatures we can use the non rel¬ 
ativistic approximation so that our problem is equivalent 
to the classical computation of the viscosity of a gas of 
hard spheres of radius R so that a = By using the 
well known result in this case (see for instance Q) we 


find the viscosity to be given by 
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This formula provides a nice estimation of the viscosity 
of the pion gas. However it does not take into account 
any quantum effect that can be important at low temper¬ 
atures. In addition it does not include any dependence 
in the pion chemical potential, i.e. in the pion density. 

In the following we will show how these effects can be 
properly taken into account to find the pion gas viscos¬ 
ity in terms of the temperature and chemical potential 
(or density) in the low temperature and density regime 
at the leading order. We will also estimate the mag¬ 
nitude of the next to leading corrections. The plan of 
the paper goes as follows: In sec.2 we review the main 
equations of the non-relativistic ideal pion gas to set the 
notation and for further reference. In sec.3 we describe 
briefly the fundamentals of the kinetic theory needed to 
follow the present work and in sec.4 we show how the 
hydrodynamic equations emerge from it. In sec.5 we ob¬ 
tain the equations for the transport coefficients in terms 
of the appropriate solutions of the linearized transport 
equations. In sec .6 we simplify these equations to write 
them in useful form for finding explicit formulae for the 
viscosity. In sec.7 we give the cross-section for the low 
energy elastic pion scattering obtained from the lowest 
order chiral perturbation theory. In sec .8 we obtain the 
lowest order terms of the pion gas viscosity. In sec.9 we 
show our numerical results and discuss the magnitude of 
the next to leading corrections and in sec. 10 we set the 
main conclusions of our work. Finally in the appendix 
we study the properties of the polynomials introduced in 
this work that play a role similar to the standard Sonine 
polynomials in the more familiar classic computations. 


II. STATE EQUATION 


In this section we review briefly the equation of state 
of a free non-relativistic bosonic gas, which will provide 
our description of the pion gas at low temperature and 
density. As it is well known [Q, for a large volume the 
pressure can be written as 


P= - 


9'^ r .1 2 , 

^ log [1 - eSo] 


1 _ g/3(A‘-£^(p)) 


( 3 ) 


where g is the number of pions (5 = 3), T = 1//3 is the 
temperature, V is the volume, /i is the pion chemical po¬ 
tential and E{p) is the energy of a non relativistic pions in 
terms of the momentum i.e. p = E{p) = p^/2M.n.. 

The number density or the number of pions per unit of 
volume is given by 




dp p^ 


I 

e-0ii^-E{p}) — 1 


9 1 

Ve-dM - I 


( 4 ) 
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where the last term is the number of pions with zero 
momentum per unit of volume. When this term is not 
negligible it is responsible for the Bose-Einstein conden¬ 
sate. We thus define 


_ g 1 
V e-^f^ - 1 


(5) 


where the integrals can be written in terms of the gn 
functions as 



dx = r(n -I- l)gn+i{z) 


( 10 ) 


with 


Clearly the above equations make sense only when the 
chemical potential satisfies fj. < 0. As a consequence, in 
the thermodynamic limit, where N^V —>■ oo with N/V 
constant, we find two phases: In the first one (normal 
phase) fJ. < 0 and hq = 0. However if we lower the tem¬ 
perature, keeping a fixed density, g increases until g = 0. 
At that point no 7 ^ 0 and the ground state density starts 
to grow forming the Bose-Einstein condensate. Eventu¬ 
ally, at T = 0 , all pions are in the ground state so that 
n = ng. The critical temperature Tc where the phase 
transition occurs can be found to be 

' V9C3/27 

where (' 3/2 is the Riemann C function evaluated at 3/2. 
Therefore the critical temperature can be made arbitrary 
small by lowering the density. In this work we are inter¬ 
ested in the the low temperature and the low density 
regime. In order to simplify the analysis we will not con¬ 
sider the contribution of the condensate. Thus for a given 
temperature T, we will always assume that the density 
will be low enough so that T > and no significant 
fraction of pions will have zero momentum (see for a 
discussion of the pion condensate). 

In the following it will also be convenient to introduce 
the fugacity 


z = exp 


(f) 


and then the above equations for the density and the 
pressure can be written as 


g 

27r^ 


P 


'0 — 1 


-dp 


( 6 ) 


and 


P = - 


27r^ 


In ( 1 — ze 


dp (7) 


Introducing the variable x = p^/{2 MtjT) we have 




( 8 ) 




1=1 


ln+1 


SO that 


gs(l) = Cs 


and 


gs-i{.z) = z 


dgsjz) 

dz 


Therefore we have for a non relativistic free pion gas at 
equilibrium at low temperature and density 


and 


n = g 


P = g 


M^T 

27T 


m^t 

27T 


3/2 


g3/2 (z) 


3/2 


Tg5/2 {z) 


Thus the equation of state can be written as 

.g5/2 {z) 


P = nT 


g3/2 {z) 


( 11 ) 


( 12 ) 


(13) 


III. THE KINETIC THEORY FOR PIONS 


The possibility of having an hydrodynamic description 
of the pion gas requires the definition of macroscopic 
fields such as the velocity V = V {r,t), the pressure 
P — P (r,t), the number density n = n{r,t), the temper¬ 
ature T = T {r, t) and others. In other words, it is needed 
to have a large number of pions in each volume element 
dV. The statistical description of the gas is based on the 
(one-particle) distribution function / (r, n, t) which gives 
the number of pions inside the phase-space volume ele¬ 
ment dr dv at the instant t. In fact the phase-space vol¬ 
ume element is given by dr dp but in the case considered 
here the difference amounts to a constant factor which 
we include in the definition of the distribution function 
/. Thus the total number of pions will be given by 

N = J f {f,v,t) dfdv = j n{f,t)dr (14) 


and 

P=£^ T ^ZT^P^dx (9) 


As we are interested here in the behaviour of the pion 
gas at low temperatures and densities, most of the the 
pion interactions will occur at low energies. Then the 
pion-pion interactions will be essentially elastic and the 
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total number of pions will be a conserved quantity. This 
in particular means that it makes sense to introduce the 
chemical potential /x = /i (r, t) related with the pion num¬ 
ber even in presence of low-energy pion interactions. The 
evolution of the distribution function is determined by 
the BBKGY hierarchy (from N.N. Bogoliubov, M. Born, 
G. Kirkwood, H.S. Green and J. Yvon) l^. This is a set of 
N coupled equations for the n-particle distribution func¬ 
tions /(„) {fi,vi,f2,V2,—,rn,Vn,t) with / = /(i). How¬ 
ever, for low-density gases it can be a good approxima¬ 
tion to truncate the BBKGY hierarchy. For example by 
using the Bogoliubov ansatz, i.e. by writing the two- 
particle function in terms of an appropriate product of 
two one-particle functions (molecular chaos hypothesis), 
it is possible to decouple the first BBKGY equation from 
the others to obtain the well known Boltzmann equation. 
This equation was modified by Uehling and Uhlenbeck in 
order to include the quantum statistic effects (see for ex¬ 
ample l]^). For the bosonic case the Uehling-Uhlenbeck 
equation reads 


(a* + FiV,-) h = c [/i] 


where 

C[fi] 


= / dv2dau 





-hh (i + 4) (i +1) 


(15) 


(16) 


is the collision term or collision functional. The above 
equation describes the irreversible space-time evolution 
of the (one-point) distribution function /, being this ir¬ 
reversible behaviour a consequence of the truncation of 
the BBKGY hierarchy. As usual we use for short the no¬ 
tation fi = f {ri,Vi, t) and // = / (rj, F', t) with i = 1,2. 
da is the differential cross-section corresponding to the 
elastic scattering of two pions with initial and final veloc¬ 
ities Fi, V 2 and F^, F 2 , respectively, u is the relative veloc¬ 
ity of the initial particles m =| Fi — F 2 | and, finally, the 
normalization constant A is given by A = gM^/ (27r)^. 

As it is the case of the Boltzmann equation, the 
Uehling-Uhlenbeck equation typically drives the gas from 
arbitrary initial distributions to the equilibrium distribu¬ 
tion which is a fixed point of this equation in the sense 


J 


that the collision functional vanishes on it 

C[f]=0 


(17) 


For a gas moving at a macroscopic velocity V the equi¬ 
librium distribution corresponding to temperature T and 
pion chemical potential g is given by 


f (F, v,t)= Ale 




-1 


- 1 


(18) 


from which the properties of free pion gas given in the 
previous section can be derived. The case of local equi¬ 
librium is also described by this distribution but having 
the V, T and /i parameters depending on the position 
and time. This case corresponds to an ideal fluid. How¬ 
ever, in order to describe dissipative processes, such as 
viscosity or thermal conduction, it is needed to consider 
the case where the pion gas is not even in local thermal 
equilibrium. 


IV. HYDRODYNAMIC EQUATIONS 


Given some magnitude ip it is possible to compute its 
statistic average by using the distribution function as 


/^\ = I ^ I fi’dr 

f fdf n (r, t) 


(19) 


where n = n{r, t) is the pion number density. This aver¬ 
age can be understood as the corresponding macroscopic 
magnitude in the hydrodynamic approach |pd| . For ex¬ 
ample, the macroscopic fluid velocity V is given by 


U ^ I fvdv ^ 1 

f fdv n 


fvdv 


( 20 ) 


The equation governing the evolution of any of these 
magnitudes (transport equations) can be obtained by 
multiplying ipi by the Uehling-Uhlenbeck equation and 
then integrating Fi so that 


J dvi (dt + viVfitpi = J dvidv2dautpi {f[f2 {1 + A ^fi){l + A ^ f2) - fif2 + A ^f[){l + A V 2 )) (21) 

The r.h.s. term can be written in a more symmetric fashion as 

i J difidv 2 dv 3 dv 4 I T |2 (V^i + V^2 - V'i - ^2) (/(/2 (1 + FUVi) (1 + FI-V 2 ) - /1/2 (1 + A-VO (1 + A-VO) (22) 

where | T p is defined through the relation _ 


I T P 

da = -0 - 

\VI-V2\ 


( 23 ) 


From this form of the collision terms it becomes apparent 
than it vanishes for any quantity tp which is conserved 
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in the collision such as the total energy or momentum. to write the velocity as IT = y + ^, i.e. ^ is the pion veloc- 

On the other hand, by using the definition of statistical ity relative to the macrocopic velocity V and therefore it 

average, the first term of ( |^ ) can be written as: averages to zero. Choosing tp = l^v and (correspond¬ 

ing to the conservation of a constant, the momentum and 
dtn {ip) +Vpn {ipv) (24) the energy in the elastic collision) we obtain the equations 

Thus for conserved magnitudes this term vanishes. In or¬ 
der to find the basic hydrodynamic equations it is useful 

I 


ip = 1 —> dtn + VnV = 0 

tP = V —> dtnV^ + + V,n (^^j) = 0 (25) 

tP = ^ dtn {V^ + (e")) + {V^ + (f)) + V,2nC, (^.0) + V.n ) = 0 


In order to understand the meaning of the above equa¬ 
tions we can introduce the following macroscopic quan¬ 
tities: 

P„-=M,n(e.C,) (26) 

is the pressure tensor, 


In spite of the above equations for the ideal fluid, it 
is well established that an irreversible flux of energy and 
momentum appear in an inhomogeneous gas. This flux 
gives rise to the well known dissipative processes of ther¬ 
mal and conduction and viscosity. For example an small 
gradient of temperature produces a heat flow which in a 
first approximation can be written as 


Q = lM^n{e) (27) 

is the internal energy density and 

Qz = (28) 

is the energy (heat) flux. Due to the isotropy of the pion 
fluid the pressure tensor can be taken as diagonal, i.e. 
Pij = P6ij. In this case Q = 3P/2 and qt = 0. Thus, (by 
taking V = 0 but not its derivatives, which can always 
be done) we arrive to the standard ideal fluid equations, 
namely the continuity equation 


dtn + n\/iVi = 0, (29) 


the Euler equation 


= -K^ipr (32) 

where k is a coefficient called thermal conductivity. In 
order to include momentum diffusion effects we add an 
extra term to the isotropic pressure tensor: 


P.,=PS.,+Pl^ (33) 


The requirement for this tensor to vanish for uniform 
translations and rotations of the fluid, makes it possible 
to write it to first order in the velocity derivatives as 


l dxj dxi 



cE 


dV, 

dxi^^^ 


(34) 


nM^dtV, -H = 0 (30) 

and the energy conservation equation 

dtQ + lQV^V, = 0 (31) 


Thus it is needed to introduce two coefficients rj and ( 
which are usually called first and second viscosity coef¬ 
ficients. The introduction of the dissipative effects does 
not modify the continuity equation. However the Euler 
equation becomes 


J 


dtnVi -I- VjnViVj + -j^^jPSij 



V 


dxj dxi 



+ cE 


dxi 


(35) 
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which is known as the Navier-Stokes equation. The energy conservation equation is in this case 


dt ( nV^ + J-Q ) + V.F. ( nV^ + = 


ViKViT + T^ViV,- 


dVi dVj z A 

fi'T ■ r)'r: ^ r)'ri 


2 

M, 

3 


2 

"Mi 


dVi 


dxj dxi 


3 dxi 
1=1 


9V; 


/ = 1 


(36) 


Thus the description of the dissipative flow of energy and 
momentum requires the introduction of three transport 
coefficients k, rj and (^. However, for an incompressible 
fluid the second (shear) viscosity does not play any role. 
In the following we will assume that this is the case for 
the pion gas considered here. 

V. COMPUTATION OF THE TRANSPORT 
COEFFICIENTS 

As it has been mentioned above the description of the 
dissipative processes requires going beyond local equi¬ 
librium. Thus, in order to compute the transport co¬ 
efficients we have to consider a distribution function / 
slightly perturbed from the equilibrium distribution /°, 
i.e. 

f = f + Sf = f + f^ (37) 

where x is an arbitrary function of the velocities which 
represents the inhomogeneous contribution to the distri¬ 
bution function. This contribution will be assumed to be 
small in the sense that x/T ^ 1) since we are interested 
only in the computation of the transport coefficients. The 
X function must by determined by solving the Uehling- 
Uhlenbeck equation which can be linearized with respect 
to X- Then it is possible to compute the transport coeffi¬ 
cients from this function. For example, the heat flux can 
be written as 

Qi = t) 


For an isotropic gas Kij = nSij so that, we have the fol¬ 
lowing equation for the heat conductivity 

« = ^ / fhv'^dv (43) 

In a similar way it is possible to write the viscosity in 
terms of an h function. In order to do that we write the 
pressure tensor in terms of the distribution function 

Pij = J t) 

For the present calculation it is enough to consider the 
case V = 0 (but not its derivatives). Then the non 
isotropic component of the pressure tensor can be written 
as 

= ^ / fxViVjdv (45) 

On the other hand for an incompressible fluid this tensor 
is 

Plj = -2r? (46) 

where we have introduced Vij = {diVj +djVi)/2. For 
further convenience we choose the perturbation to have 
the form 

X = - (Vki - hki {v) (47) 


Taking into account that this flux vanishes for an homo¬ 
geneous gas and using the definition of the heat conduc¬ 
tivity we have 

-kV*T = ^ y fxviv^dv (39) 

Now it is useful to write the x function as 


where hki is a microscopic velocity dependent quantity 
which can be written as 


hki 


VkVi - -v'^Ski ) h {v) 


(48) 


Therefore, by equating the two forms of P/- we find 


X = -'^iThi {v) = -\/iTvih (v) 

Therefore 


(40) 


2 ? 71 Vij - -VhVhS^j 


'Hij kl 



1 


^ hVhdki 


(49) 


kV iT = KijV jT 

with the tensor conductivity Kij being given by 
K^j = ^ f fhvjViv'^dv 


(41) where ijijki (the viscosity tensor) is given by 


(42) 


rjijki = M^T J VtVj (^VkVi - /° (v) h (v) dv 


(50) 



7 


Due to the several symmetries of this tensor it can written and the distribution function as 
as 


Tlijkl — Tj 





(51) 


so that we have 

V = ^ J v^VJ (^v^Vj - f (n) h (v) dv (52) 


f 

T 


/ = /° + 


Uj 




s =0 




(57) 


As the h function depends only on the velocity modulus 
it is possible to perform the angular integrations 


V = 


47r 


A 


z 1 exp 


(^-i) 


v^h{v)dv (53) 


In the following it will be useful to consider h as a func¬ 
tion of the adimensional variable x, i.e. h = h {x) defined 
as 


For all the families of polynomials, i.e., for all the r val¬ 
ues, we have exact conservation of the collision invariants. 
This fact makes possible to choose r so that we get the 
simplest integral. The appropriate value turns out to be 
r = ^12 and then we get 


_ 

“ 2T 


(54) 


At this point it is customary to develop the h (x) in terms 
of the Sonine polynomials. However it is more appro¬ 
priate in our case to introduce a new family of orthog¬ 
onal polynomials (z; x) since we are considering the 
Uehling-Uhlenbeck as the transport equation (i.e. we 
have taken into account quantum effects) instead of the 
more common Boltzmann equation. The definition and 
the main properties of this new polynomials can be found 
in appendix A. In terms of the h function can be writ¬ 
ten as 


A7r3/2 


2,T \ 

JTj 97/2{z)B, 


(58) 


In this way the computation of the viscosity of the pion 
gas has been reduced to the computation of the Bq coef¬ 
ficient. 

VI. SIMPLIFYING THE 
UEHLING-UHLENBECK EQUATION 


h{x) =^BsPr (z]x) (55) 

s=0 


Therefore 


lOT 



OO 


Y.b,p:{z-,x) 

(56) 


As discussed above the the perturbed distribution 
function must be a solution of the Uehling-Uhlenbeck 
equation (ITsI). In order to solve this equation it is quite 
convenient to write the l.h.s. in terms of macroscopic 
quantities, as the temperature, velocity, density and fu- 
gacity, and their derivatives. This can be done by using 
the equation of state and the hydrodynamic equations as 
follows: First of all at lowest order we have 


J 


dtf + v^f-dtf + dVf = Adt 

Thus the temporal derivative is 


AvV 


z-iexp(f)-l/ lz-iexp(f)-l 


/°Viexp(f) 


dtf = - 


A 


E 


-dtE-—dtT-dt\ogz 


where the energy E is given by 




") 


2 


(59) 


(60) 


(61) 











so that we have 


dtf° = 




In a similar way we can obtain 




Azexp (—y) 


^ ^ V, log^ + ^VjT + 


T 


Thus we have 


dtf + VjVjf = 


f 


Az exp (— 


fr. E ^ ^ E \ 

logz + jT^dtT + —v^^tVi + VjWj logz+ j^VjVjT + —VrVjVijj 


(62) 


(63) 


(64) 


Here it is possible to use the state equation for the free pion gas and the ideal fluid equations (Euler, continuity and 
energy conservation) to find 


dtf + 


r [1 

Hz exp (—y) T 



5 55/2_W\ 

2 53/2 (.z) J 


ViViT + 



ViVj 




(65) 


On the other hand, substituting expression ( |3^ into expanding to first order in the perturbation y and taking 
benefit from the equilibrium distribution function relation 


/07r(l + H-Vi°) (l + H-i 


/2°) =/?/2° (l +/?') (l + H-i/o') 


( 66 ) 


we obtain 


/(/' (1 + H-i/i) (1 + H-V 2 ) - / 1/2 (1 + H-i/O (1 + H-V 2 ) = 


H 2 


.-^pE/T 


T - 1) {z-^e^^/^ - l) {z-^e^'E^ - l) - l) 

where E = Ei + E 2 = E[ + E 2 and the A symbol is defined as 

^ [/ (a^)] = / (a;i) + / {X 2 ) - f (xi) - / (a;2) 

Thus the whole collision term can be written as 


X ( 1 - zexp ( -p 


C [/i] = J dv 2 da uz ^ exp {E/T) /i/ 2 /(/ 2 A 


X 1 - z exp - 


(67) 


( 68 ) 


(69) 


where we have omitted the superindex 0 in the equilibrium distribution function. In this way the complete transport 
equation for the slightly unhomogeneous pion gas becomes 


n 


ze 2 t 


Mt,vI _ 5 55/2 jz) 
‘ET 2 33/2 {z) 


J dv 2 da u z +“ ) f^f 2 f[f 2 A 



Vij 


M„v' 

1 — ze 

2T 


(70) 


where we have defined U = {vi + V 2 ) /2 = {v[ + V 2 ) /2, 
i.e. it is the collision center of mass velocity. Once again 
M = ui — U 2 is the relative velocity of the incident pions. 
Thus, in particular, we have 



217" + -E 


(71) 


VII. LOW-ENERGY PION CROSS-SECTION 


In order to solve the above transport equation it is 
needed to know the differential cross section for pion scat¬ 
tering. As we are interested in the low energy regime 
the most appropriate approach is the Chiral Perturba- 
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tion Theory which relies in the chiral symmetry pattern 
of the strong interactions. In particular it is quite useful 
the chiral Lagrangian approach suggested by Weinberg 
and extended up to the one loop level by Gasser and 
Leutwyler. This chiral Lagrangian approach consists in 


J 


a systematic expansion on the pion field derivatives and 
the pion mass which are considered to be of the same 
order. From the lowest order of this expansion, the pion 
elastic scattering amplitude is given by 


■-Ml 


t- Ml 


•-Ml 


^ai,a2,aya2 p2 ^aia2^a[a2 “f p2 ^aia[^a2a2 Jp2 ^ 010 ^^ 020 ^ 

-^TT -^TT -^TT 

I- 


(72) 


where the standard Mandelstan variables s, t and u are 
related hy s + t + u = iMl, the a subindices refer to the 
pion isospin and is the pion decay constant. In the 
total isospin basis there are three independent amplitudes 
corresponding to I = 0,1, 2 which are given by 


^0 _ 


= 

rp2 _ 


2s-Ml 
F2 

TT 

t — u 


^2 
-*■ TT 

2M1-S 

F2 


(73) 


The averaged modulo squared amplitude is given by 


\Tf = 


which leads to 


E/ (2/+1) 


^(2/+l)|T^ 


(74) 


17^1' = ^ {mmI + 9s^ - 2AMls + 3{t- uf) (75) 

The relativistic differential cross-section for elastic pion 
scattering is 


da = 


\T\^ 


4:uEiE2 


-dLIPS 


(76) 


where the differential Lorentz Invariant Phase Space is 
given by 


dLIPS = C {2Trf (pi +P2- p'l - P 2 ) 


dp'i 


dp '2 


2 (27r)'^ E[ 2 (27r)^ S' 


(77) 


where the constant (7=1/2 takes into account the iden¬ 
tity of the final pions. 

As discussed above we are interested in the non rel¬ 
ativistic regime where pi = M^^Vi, P = 2M,rC7 and 
Ei = Mtt + MTrvl/2. In this limit it is convenient to 
write the cross section in terms of the center of mass and 
relative velocities U, u and u' defined as 


(78) 


After some standard calculations the non relativistic re¬ 
duction of the cross section can be found to be 


^ 23 Ml / 3^ . 13 

da = -j ( 1- ^ - u 

768 TrSi V 2 184 


2 rr2 „,2U + ucos6 \ , 

^ -f UMosO' --- dcosO' 

2U cos 9' + u J 


where 9' is the angle between by U and u'. 


(79) 


VIII. LOWEST ORDER VISCOSITY FOR THE 
PION GAS 


As we have seen in previous sections the perturbation 
of the distribution function appropriate for the compu¬ 


tation of the pion gas viscosity can be written as 


1 . 


1 


X = -[Vij - - -V 4 

00 

Y,B,Pi,^{z-x) (80) 


s=0 
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where the Bg coefficients must be obtained by solving the 
linerized transport equation being Bq the dominant one. 
In order to do so let us consider for a while a tensor l^i 
which is an arbitrary function on x and Di. For further 
convenience we have defined the hat on any velocity as 


V = 


2T 


(81) 


By multiplying this tensor by the transport equation and 
integrating on vi we have 


dvilki iz;vi) 


f? 


vl- 




1 /* 

/ dvidv2da lu {z\Vi)u 


(l - ze (z; v^) 


(82) 


r 


where we have assumed the divergence of the velocity to In order to compute the Bg coefficients we have to pro¬ 
vanish since the pion gas is taken to be incompressible. ceed in a similar way with the r.h.s. of the equation (|82|). 


Now we choose the tensor to be 


Iki {z;v) = (^VkVi - ^v'^Ski^ (l- ze " ^ P 5/2 (■ 2 ;^^) 

(83) 

Then, one of the two integrals appearing in the r.h.s. of 
the equation (§2|) vanishes, namely 


For that purpose we will take the low energy limit of the 
differential cross section which is just 


7687rFi 


( 86 ) 


dvilki {z-,vi) 


n 


ze 1 


5 95/2 {z) 
2 53/2 (z) 


wiiViT = 0 


(84) 


where 6' belongs to the interval (0, tt). Therefore, at very 
low energies the elastic pion scattering cross section can 
be considered as a constant, so that it is formally equiv¬ 
alent to the elastic scattering of rigid spheres of radius 


Therefore, the two sides of the equation (^ are pro¬ 
portional to Vij which is arbitrary. Thus one possible 
solution to the equation is found by eliminating this ten- 


R = 


/ 23 


ttF3 V 384 


(87) 


After the manipulations described below the r.h.s of the 
sor in both sides and then by contracting with the tensors transport equation can be written the collision term as 


dikdji so that the r.h.s. of the equation (p2D becomes 






231^ 

768F4 \M^) ^ ^ 


s=0 


/ 2T \ 

Ur 57/2 (.) ^ot (85) where we have defined 
Z \ IVi-jY J 


( 88 ) 


/*00 nlT n'2'K 

hst{z) = / dUdu / dO'dO / (i(/>sin^'sin^ TT ( 

Jo Jo Jo 

VkVi - (^1 - ze"'"'') P 5/2 (z; 


z-ie< - 1 


-1 


dl ,a2 

VkVl - ^ 4 ;!’^) (1 - 26"" ) P 5/2 {z]v'^) 


(89) 


Then the transport equation (p2D can be written as 


Therefore by defining 


23M^A / 2T 


19207r3/2F4 


97/2 (^) = '^ 0 * 


s=0 


23AUA 

= 19207r3/2p4 ) 57/2 (^) Bs (91) 
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the viscosity is given by 


IX. RESULTS AND DISCUSSION 


77 = 


2T 9607r3F4 ^ 


Mt, 23M^ 

where Bq is a solution of the equation system 


'^Bsbts (z) = 6to 


(92) 


(93) 


s=0 


or 


fbooiz) 601 ( 2 ) 602 ( 2 ) 
610 ( 2 ) 611 (z) 612 ( 2 ) 
620 ( 2 ) 621 ( 2 ) 622 ( 2 ) 

V : ; 


\ 

fBo\ 




Bi 


0 


B 2 


0 

/ 



W 


(94) 


In order to find the leading contribution to the viscosity 
we can truncate the above system to find 


Bo = 


1 


boo ( 2 ) 


(95) 


In this way the calculation of the viscosity have been 
reduced to the computation of the integral 690 ( 2 ). 


According to our previous discussion the viscosity of 
the pion gas at low density and temperatures is given by 
equation (|^), where Bq is the solution of the equation 
system given i n (|93| ) and the corresponding coefficients 
are defined in (|89[). Thus it is possible, at least in princi¬ 
ple, to find the viscosity as a function of the temperature 
T and the fugacity z (or the pion chemical potential). 
Alternatively it is also possible to find the viscosity as a 
function of the temperature and the pion number den¬ 
sity. Due to the complicated integrals appearing in the 
definition of the bgt ( 2 ) functions it has been possible to 
get numerical results only by Monte Carlo |Q integra¬ 
tion. For the sake of simplicity most of these numerical 
results have been obtained by using the approximation of 
( P^ ) and the hard sphere cross section of ( 86 ). In order 
to check the validity of these approximations we have 
computed the 4x4 matrix containing the first bst ( 2 ) 
functions for the particular value z = 1. The result and 
the corresponding errors is 


bst. — 


/ 277.2 ±0.7 I -42 ±1 -99 ± 3 -408 ± 14 \ 

-42 ± 1 1043 ± 8 -395 ± 8 -800 ± 30 

-99 ± 3 -395 ± 8 10700 ± 110 -7100 ± 200 

-408 ± 14 -800 ± 30 -7100 ±200 199000 ± 2000 / 


(96) 


Now we consider the submatrices 1x1,2x2,3x3 and 4x4 and solve the corresponding system for the appropriate 
number of B so that we can check the convergence of the results as the size of the truncation is increased. The results 
are shown in table 1. There we see that the error is of the order of 1.6 % 



Bo 

Bi 

B 2 

Bo 

Zero order 

First order 

Second order 

Third order 

3.608-10-3 

3.630-10-3 

3.647-10-3 

3.666-10-3 

1.497- 10-4 
1.654- 10-4 
1.765-10-4 

3.985 - 10-3 
4.700 - 10-3 

9.902 -10-3 


To take into account the velocity effects in the cross section we will consider (79) instead of ( p 6 [ ) in the computation 
of the bij integrals. Then we will have 




where 6 °- is the leading term corresponding to the hard sphere approximation and 


^00 — 


dUdu / dOda 


d4> 


sin 6 sin a + 2 “ 


13 , 2U + ucos9 

—V? - 31/2 21/2 cos 9— --- 

92 21/ cos 9 + u 


0 - 1 ) - 1 ) (e"? - l) - l) 

1 


0 ( 


1 — e 


(97) 


(98) 
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At the zero order in the sense of the previous table we 
have 


Bn = ^ = 


boo 6E]o + ^Ko 
which at low temperatures can be written as 




so that 

V = 


2r 9607r3F4^2 1 A T 

97/2 (^) ^0 I 1 ^50 

Coo \ "00 


23M, 

This formula can be applied provided 


(99) 


( 100 ) 


( 101 ) 





( 102 ) 


After a numerical integration we find, for z — 1) ^00 ~ 
—310±40. By comparison with our previous result 5 qq = 
277.2 ± 0.7 we get the upper bound for the temperature 


T< (0.90 ±0.12) (103) 


which means that the hard sphere approximation is safe 
provided the temperature is much smaller than the 90% 
of the pion mass. As we are working in a non relativis¬ 
tic framework this will always be the case and then no 
additional constrain on the applicability of our results is 
coming from the hard sphere approximation. 

Therefore in the following we will work at zero order 
and in the hard sphere approximation. In this case the 
viscosity is given by 


V (T, z) 


^ 4096 ^ 2 ^ 

345M^/^ 


{boo {z))-^ 



z~^e^ — 1 



(104) 


Some numerical results obtained from this formula are 
found in table 2 for different values of the fugacity z. 



g (MeV^) 

Quantum case (z = 1) 

^/^4.00• 10® 

Quantum case (z = 0.05) 

Vt 1.40 • 10® 

Classical case 

^/T 1.9- 10® 


The classical case is just the limit of z going to zero. 
This corresponds to a classical gas of hard spheres of 
radious 


fW 
ttF^ V 384 


(105) 



0 0.2 0.4 0.6 0.8 1 


FIG. 1: Plot of (&00 (z)) ^ versus z. It can 

be understood as a plot of the viscosity, for constant temper¬ 
ature and in arbitrary units, versus the fugacity. 


Because of the particular way we have performed our 
computations the particular case z = 0 is not numerically 
accessible. However this case has been treated long time 
ago and it is possible to find an analytic result for the 
viscosity, namely 


77 = 


/— 

167^772 


(106) 


Thus we can check that our computations have the 
proper z = 0 limit. In fig.l we show the behaviour of boo 
and the viscosity, for constant T, in terms of the fugacity. 
From this plots we learn that the viscosity of the bosonic 
quantum gas is smaller than that of the classical gas. 
This could be expected because of the following heuristic 
argument: After some elastic collision the emerging par¬ 
ticles have more affinity for occupied states than classical 
ones. Therefore the microscopic transfer of momentum 
is more effective in a classical gas than in a bosonic gas 
and then the viscosity is also larger. A pictorial view of 
this fact is displayed in fig.2. 

In many cases it is useful to have the viscosity as a 
function of the temperature and the density instead of 
the temperature and the fugacity. The pionic density can 
be written in terms of the temperature and the fugacity 
as 




(107) 


This equation defines also implicitly the fugacity as a 
function of the density and the temperature z = z {n,T). 
However, due to the complexities of the above integral 
this functions have only been computed numerically. In 
fig.4 we plot ri/\/T as a function of The com¬ 

puted points can be fitted quite well in the plotted range 
with a polynomial 


/i {x) = A + Bx ± Cx^ 


(108) 
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o o 
o o 

O o 
O 

o o 


with constants A = 0.1814±0.0004, B = (—2.42 ± 0.03) • 
10“^ and C = (9.6 ± 0.3) • 10“®. Thus the function ob¬ 
tained (fig.3) for the viscosity versus the temperature and 
the density is 


o o 
o o 

O o 

O 

o o 


Q 


o o 
o o 
O o 

O 

o o 


FIG. 2: a) Two particles, belonging to two different regions of 
the gas with different average momentum (black and white), 
before an elastic collision, b) Particles before the collision in 
the classical case where they typically get a new momentum 
(represented in grey) c) Particles in the bosonic gas where 
they typically want to go to occupied states with the same 
momentum. 


J 


V {T, n) 


4096 ^ 2 ^^ 

345Mi?/^ 


0.1814 


2.42- 10"'‘^^ -b9.6- 10"® 



(109) 


0.18 I-1-1-1-1-1-1-1-1 

0.16 - 

0.14 - V 

0.12 - 

0.1 - 

0.08 - 

0.06 - 

0.04 I-'-'-'-'-'-'-'- 

0 100 200 300 400 500 600 700 800 

FIG. 3: Numerical fit of rj/\/T oc /i ^ 


As it may be noticed in fig.4, the values of the viscosity 
cover the range from 10® to 10^ MeV® for the considered 
densities and temperatures. 

X. CONCLUSIONS 

In this work we have computed the viscosity of the pion 
gas, starting from first principles only. This computation 


is relevant for the hydrodynamic description of hadronic 
matter at low energies and densities. The work relies on 
the use of the Uehling-Uhlenbeck equation, which is the 
quantum version of the Boltzmann equation, and chiral 
symmetry, which determines completely the structure of 
the pion interactions at low energies. We have also found 
a formula for the viscosity in terms of the temperature 
and the density which properly fits our numerical results. 
The interest of the computations is twofold. First it is 
useful to check when the usual assumption of perfect fluid 
for the hadronic fluid is reasonable enough. Second it 
can be used to include the viscosity in the Navier-Stokes 
equations when dissipative effects cannot be neglected. 
The main limitation of our work is that it can be ap¬ 
plied only in the non-relativistic domain. As a conse¬ 
quence of that it cannot be used directly in the study 
of the events observed in the modern ultra relativistic 
heavy ion colliders such as RHIC or the future LHC. In 
any case we understand that our work is interesting in 
order to show that the computations of transport coeffi¬ 
cients of the pion gas can be done in a complete model 
independent way. In fact we consider that the result pre¬ 
sented here is just the first step of a complete relativistic 
computation which could be applied at higher energies 
and therefore in more realistic situations, to improve the 
standard hydrodynamic description of hadronic matter. 
Work is in progress in that direction. 
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FIG. 4: 3D plot of viscosity vs. temp erature and particle 
number density, given by equation (|l09|). 
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XI. APPENDIX 

In this appendix we study the main properties of the 
family of orthogonal polynomials P® (z; x) defined on the 


interval ( 0 , oo) with measure 

where z G (0 1] and r > 0. By definition the polynomials 
are orthogonal for some given r, or in other words, 

(p;, P®') = dx ^_y^"_^ P® (z; x) P: (z; x) 

= A(z;r, s)(5ss/(lll) 

For simplicity we define the polynomials so that they are 
monic, i.e. the coefficient of the term of highest degree 
in each polynomial is taken to be one. Thus, the first 
polynomial is always the unity and his norm is given by 

p; = 1 ^ (p;, p,») = 

= r (r -f 1 ) (z) ( 112 ) 

From the condition 

(P°,P;)=0 (113) 

it is possible to compute the second polynomial which 
turns out to be 

(z; x) = (r -bl) - a; (114) 

gr+l [z) 

The third polynomial is obtained from the conditions 
(P 0 ,P 2 ) = 0 and (P^\P 2 ) = 0 which gives 


Pr (^; x) 


(r P 3) gr +2 (z) gr+i (z) -{r + 2) g^^^ (z) 

(r + 2 ) g.+i (z) gr +3 (z) - (r + 1 ) 5,^2 (z) ^ ^ ^ ^ ^ 

{r + 3) gr+i (z) gr+4 {z) - {r + 1) gr +2 {z) gr+3 {z) , , ‘ 

{r + 1) gr +2 (z) - (J' + 2) gr+i (z) gr+3 (z) 


Higher polynomials can be obtained in a similar way. 
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